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ABSTRACT 


The  non-cooperative-target  beam-control  problem  has  been  the  subject  of  intense 
investigation  since  the  synthetic-beacon  concept  was  first  introduced  to  the  high-energy-laser 
community  in  1982.  While  numerous  analytical  studies  and  computer  simulations  have  been 
performed  to  evaluate  the  practical  utility  of  this  phase-measurement  technique,  prior  to  Lincoln 
Laboratory’s  SWAT  (Short-Wavelength  Adaptive  Techniques)  program  no  experimental 
verification  had  been  obtained.  In  the  first  phase  of  the  SWAT  investigation,  completed  in  1985, 
a  high  degree  of  correlation  between  differential-phase  measurements  from  natural  and  artificial 
sources  was  demonstrated. 

The  next  phase  of  the  SWAT  program  will  be  performed  at  the  AMOS  (Air  Force  Maui 
Optical  Station)  facility  in  Maui,  where  a  241 -actuator  adaptive-optics  system  and  an  array  of  six 
dye  lasers  will  be  integrated  with  the  site's  60-cm  beam  director.  Prior  to  shipment,  the 
adaptive-optics  subsystem  was  subjected  to  a  thorough  laboratory  evaluation,  which  culminated 
in  a  series  of  compensation  tests  involving  simulated  beacon  sources.  The  results  of  these 
measurements  are  in  good  agreement  with  theoretical  predictions  and  provide  strong  evidence  of 
the  efficacy  of  the  synthetic-beacon  approach. 


PREFACE 


From  1982  to  1991  Lincoln  Laboratory  was  supported  by  the  Defense  Department  to 
develop  the  technology  of  uncooperative-target  adaptive  optics.  A  principal  feature  of  this 
program  was  the  development  of  synthetic  beacons  as  the  reference  source  for  the  adaptive-optics 
wavefront  sensor.  This  technique  overcame  problems  associated  with  point-ahead  and  target 
brightness. 

In  the  late  1980’s  it  became  clear  that  synthetic-beacon  technology  could  be  of  great  value 
to  the  astronomical  community,  who  routinely  deal  with  dim  objects.  Moreover  the  basic 
concept  had  been  reinvented  by  the  astronomers  with,  apparently,  no  input  from  the  US  Defense 
community.  The  civilian  community  seemed  on  the  verge  of  embarking  on  a  major  development 
program  duplicating  work  the  US  Government  had  already  paid  for. 

It  seemed  incumbent  therefore  on  the  DoD  to  provide  to  the  astronomers  the  wealth  of 
information  that  had  been  developed  over  the  past  decade.  So,  in  May  1991  the  US  Air  Force 
released  virtually  all  of  the  development  work  and  associated  documentation.  The  immediate 
impact  was  that  Defense  workers  began  to  present  their  material  in  open  publication.  We  are 
now  moving  to  republish  a  number  of  documents  that  seem  of  greatest  import.  To  some  extent 
these  documents  are  dated  as  far  as  military  work  is  concerned,  but  they  are  likely  to  be  of  great 
value  to  the  astronomers.  At  a  minimum  they  are  of  historical  and  archival  significance. 

The  report  that  follows  is  one  of  many  documents  we  are  in  the  process  of  releasing.  The 
original  report  numbered  SWP-6  (“SWAT  Laboratory  Test  Results,”  dated  16  January  1989)  is 
now  updated.  If  you  possess  the  original,  it  should  be  destroyed  in  accordance  with  current 
document-control  procedures.  We  are  providing  a  copy  of  this  report  fo  the  Defense  Technical 
Information  Center  (DTIC)  so  that  users  may  request  it  There  are  intended  to  be  no  distribution 
limitations  on  this  document. 
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1.  INTRODUCTION 


The  non-cooperative-target  beam-control  problem  has  been  the  subject  of  intense 
investigation  since  the  concept  of  using  synthetic  beacons  as  atmospheric  probes  was  introduced 
to  the  high-energy-laser  community  in  1982.  The  proposed  solutions  typically  involve  the 
placement  of  one  or  more  focused  beams  in  the  target's  point-ahead  direction  (the  so-called  "A" 
method  approach)  or  the  projection  of  a  single  collimated  beam  that  is  laterally  sheared  ("S" 
method).  While  numerous  analytical  studies  and  computer  simulations  have  been  performed  to 
evaluate  the  practical  utility  of  uncooperative-compensation  techniques,  direct  experimental 
verification  will  require  a  major  hardware-development  effort  and  an  ambitious  measurements 
program. 

To  this  end,  Lincoln  Laboratory  has  undertaken  a  series  of  field  experiments  that  are 
collectively  referred  to  as  the  SWAT  (Short-Wavelength  Adaptive  Techniques)  program.  The 
first  phase,  in  which  the  differential  subaperture  tilt  produced  by  a  synthetic  beacon  and  that 
produced  by  a  star  were  compared  in  a  pair  of  subapertures,  was  successfully  completed  in 
FY85.  The  next  set  of  experiments  will  be  performed  at  the  AMOS  (Air  Force  Maui  Optical 
Station)  facility  in  Maui,  and  will  integrate  a  241-actuator  adaptive-optics  system  and  an  array  of 
six  dye  lasers  with  the  site's  60-cm  beam  director.7  This  instrumentation  has  been  specifically 
designed  to  evaluate  single  and  multiple-beacon  A-method  deployment  schemes,  but  provision 
has  also  been  made  to  perform  a  limited  set  of  sheared-beam  tests. 

As  of  this  writing,  all  of  the  adaptive-optics  components  and  one  of  the  six  dye  lasers  for 
the  SWAT  II  experiments  have  been  installed  at  the  Maui  site.  Prior  to  shipment,  the  adaptive- 
optics  subsystem  was  subjected  to  a  thorough  laboratory  evaluation,  which  culminated  in  a  series 
of  atmospheric-compensation  tests  involving  simulated  beacon  sources.  The  stated  objectives 
included  a  demonstration  of  multiple-beacon  phase  correction  in  a  simulated  turbulence 
environment. 


1  R.  R.  Parenti,  P.  N.  Everett,  R.  Kramer,  and  D.  A.  Page,  "SWAT  II  Overview,”  Lincoln 
Laboratory  Project  Report  SWP-3  (Revision  1)  (November  1992). 
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This  report  summarizes  the  results  of  the  SWAT  laboratory  experiments  and  compares  the 
measured  performance  with  that  predicted  by  current  theory.  To  the  best  of  our  knowledge, 
these  data  represent  the  most  nearly  definitive  experimental  verification  of  the  synthetic-beacon 
concept  obtained  to  date. 
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2.  DESCRIPTION  OF  LABORATORY  EQUIPMENT 


This  section  presents  a  brief  overview  of  the  SWAT  laboratory  hardware.  A  major  concern 
in  designing  this  equipment  was  the  accuracy  with  which  the  gross  atmospheric  effects  that  give 
rise  to  focal  anisoplanatism  could  be  represented.  A  good  simulation  of  a  typical  atmosphere  was 
achieved  by  constructing  an  optical  path  having  two  sets  of  phase  screens  appropriately 
positioned  with  respect  to  the  receiving  aperture.  Also  incorporated  were  light  sources  at  scaled 
ranges  of  10  km  and  infinity,  as  well  as  a  means  to  shutter  those  sources  to  duplicate  the  action 
of  a  range-gated  beacon  laser.  All  dynamic  sequencing  functions  and  data  recording  operations 
were  performed  under  the  direct  control  of  a  microVAX  computer  system. 

2.1  OPTICAL  LAYOUT 

In  the  baseline  SWAT  field  tests,  a  set  of  pulsed-dye  lasers  will  be  used  to  project  a  total  of 
five  synthetic  beacons,  one  placed  at  the  center  of  the  receiving  aperture  and  four  others  centered 
over  each  of  four  aperture  sections.  For  passive  scoring  tests,  a  sixth  source,  usually  a  bright 
star,  will  provide  a  means  to  characterize  the  quality  of  the  phase  correction. 

All  of  the  essential  aspects  of  beacon  geometry,  phase  sensing,  wavefront  compensation, 
and  far-field  scoring  were  incorporated  into  the  design  of  the  SWAT  laboratory.  A  simplified 
diagram  of  the  optical  layout  is  given  in  Figure  2-1.  The  six  sources  representing  the  beacons 
and  the  scoring  star  were  produced  by  a  pair  of  CW  argon-ion  lasers  operating  at  5 14.5  nm.  The 
light  from  one  of  the  lasers  was  collimated;  the  other  beam  was  split  into  five  components  that 
were  laterally  displaced  to  form  the  desired  spot  pattern  and  that  could  be  individually  controlled 
with  mechanical  shutters.  In  order  to  limit  the  path  length  required  to  simulate  the  effects  of  10- 
km  sources,  the  receiving  aperture  was  demagnified  to  a  diameter  of  1.5  cm,  thus  producing  a 
scaling  factor  (designated  by  (3)  of  40  with  respect  to  the  60-cm  aperture  at  AMOS.  Scaling  laws 
for  the  critical  simulation  parameters  are  summarized  in  the  table  below. 
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Figure  2-1.  Laboratory  layout. 


TABLE  2-1 

Scaling  Factor  Summary 


Collection  Aperture  Diameter 

1/3 

Range  Parameters 

i/32 

Phase  Screen  Strengths 

pS/3 

Using  this  geometry,  the  path  length  from  the  sources  to  the  receiver  could  be  reduced  to  just 
over  6  meters. 

Wavefront  control  in  the  laboratory  was  accomplished  just  as  it  will  be  in  the  field  -  by  a 
241 -actuator  deformable  mirror  driven  by  a  modified-Hartmann  phase  sensor.  A  digital 
reconstructor  was  used  to  convert  the  phase-gradient  outputs  from  the  wavefront  sensor  to  a  set 
of  optimal  actuator  positions  for  phase  conjugation.  All  actions  were  closely  synchronized  with 
the  timing  of  the  wavefront  sensor  and  far-field  camera,  as  described  in  the  next  section. 

2.2  EXPERIMENT  CONTROL  AND  DATA  RECORDING 

The  most  significant  differences  between  the  laboratory  and  field  hardware  configurations 
are  apparent  in  the  design  of  the  control  and  recording  subsystems.  In  order  to  minimize 
complexity  and  provide  replacements  for  components  still  under  construction,  all  timing 
functions  for  the  laboratory  tests  were  handled  by  a  hard-wired  sequencer  directed  by  a 
microVAX  computer.  As  shown  in  Figure  2-2,  the  sequencer  sends  control  pulses  to  the  beam 
shutters,  tilt  mirror,  wavefront  sensor,  deformable  mirror,  and  CCD  camera  in  the  proper  order. 
The  timing  diagram  for  this  process  is  shown  in  Figure  2-3.  The  phase-correction  cycle 
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Figure  2-2.  Laboratory  experiment  control  system. 


Precursor  Shutter 


re  2-3.  Multiple-source  laboratory  experiment  baseline  sequence. 


is  completed  in  approximately  1 .2  seconds.  This  cycle  can  be  repeated  indefinitely  for  CW  loop- 
closure  experiments.  Since  the  sequencer  is  equipped  with  discrete  parallel  outputs,  die  order  in 
which  the  various  beams  are  shuttered  is  easily  modified  by  rearranging  the  output  connections. 

The  microVAX  system  shown  in  Figure  2-2  serves  both  as  a  trigger  for  the  sequencer  unit 
and  as  a  data-storage  medium.  Upon  command  from  the  system  operator,  a  data  session 
consisting  of  a  series  of  sequential  experiments  is  initiated,  with  all  essential  data  recorded  onto 
a  hard  disk.  Usually  40  experiments,  each  following  the  timing  sequence  shown  in  Figure  2-3, 
were  performed.  The  associated  records  include  the  following  information: 

•  x  and  y  phase  gradients  for  each  beacon  sample 

•  x  and  y  intensities  for  each  beacon  sample 

•  reconstructed  phases 

•  phase  drive  voltages 

•  two  sets  of  CCD  camera  images 

All  of  this  information  is  archived  according  to  date  and  time  and  is  also  immediately  available  for 
analysis  through  the  local  VAX  network.  Derived  parameters,  such  as  those  listed  below,  can  be 
readily  generated  and  presented  in  either  text  or  graphics  format: 

•  minimum,  maximum,  and  average  gradient  values 

•  gradient  curls 

•  minimum,  maximum,  and  average  intensities 

•  minimum,  maximum,  and  average  phase  values 

•  minimum,  maximum,  and  average  CCD  camera  pixel  values 

•  far-field  Strehl  ratios 
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2.3  PHASE  SCREEN  CHARACTERIZATION 

To  correlate  the  observed  performance  with  that  predicted  by  our  computer  simulations,  a 
good  characterization  of  the  phase  screens  used  in  the  laboratory  tests  is  essential.  During  most 
of  the  tests,  one  or  more  screens  from  a  set  of  three  were  placed  in  the  beam  path  at  the  positions 
indicated  in  Figure  2-1.  Digitized  phase  maps  for  each  of  these  screens  were  obtained  using  a 
high-resolution  (100  x  100)  Zygo  interferometer.  Structure  functions  derived  from  these  images, 
scaled  to  a  60-cm  aperture,  are  presented  in  Figure  2-4.  The  anticipated  shape  of  the  tilt-removed 
structure  function  for  Kolmogorov  turbulence  is  given  by  the  expression 

D y(r)  =  5 1 .3  Of2  (r/D)5/3  [  1  -(r/D)  1/3]  .  (2- 1 ) 

A  plot  of  this  function  for  a  figure  variance  of-  =  3.7  rad2  is  included  in  Figure  2-4.  All  three  of 
the  measured  curves  are  seen  to  be  close  to  the  correct  form,  with  the  match  for  Screen  #1  being 
particularly  good.  Coherence  diameter,  isoplanatic  angle,  and  log  amplitude  fluctuation  values 


for  these  screens  are  derived  from  the 

figure  variance  as  follows: 

Coherence  Diameter  = 

r0  =  0.299  D/(of2)3/5  , 

(2-2) 

Isoplanatic  Angle  = 

fr0  =  0.0942  D/h0(af2)3/5  ,  and 

(2-3) 

Log  Amplitude  Fluctuation  = 

ax2  =  9.90  (ho/kD2)5/6  0(2  , 

(2-4) 

where  ho  represents  the  altitude  of  the  phase  screen.  On  this  basis,  the  following  tabulation  is 
given.  (In  the  bench  simulation  the  screens  were  placed  in  one  of  two  positions  -  a  so-called 
low-altitude  position  that  was  adjacent  to  the  entrance  aperture,  or  a  high-altitude  position  for 
which  h0  =  2.8  km.  The  d0  and  a^2  numbers  apply  to  the  placement  of  screens  at  the  high- 
altitude  position.) 
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screen  structure  function  analysis. 


TABLE  2-2 


Phase  Screen  Parameter  Summary 


Of  (rad2) 

r0(cm) 

flo(prad) 

2  •> 
(Nepers^) 

Screen  #1 

3.7 

8.2 

9.2 

0.08 

Screen  #2 

0.8 

20.5 

23.0 

0.02 

Screen  #3 

5.6 

6.4 

7.2 

0.12 

Screens  #1  &  2 

4.5 

7.3 

8.2 

0.10 

Screens  #1,  2,  &  3 

10.1 

4.5 

5.0 

0.22 

It  can  be  seen  that  a  wide  range  of  turbulence  conditions  could  be  readily  simulated  by  adjusting 
the  number  and  position  of  the  included  screens. 
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3.  RECENT  THEORETICAL  RESULTS 


Central  to  the  SWAT  experiments  is  the  question  of  achievable  Strehl  as  a  function  of  such 
parameters  as  beacon  size,  number,  and  altitude.  Several  numerical  studies  relevant  to  the 
SWAT  configuration  have  been  conducted  in  the  last  few  years,  with  beacon-size  effects  having 
been  recently  added.  These  investigations  typically  incorporate  atmospheric  models  consisting  of 
ten  or  more  phase  screens  and  apply  geometrical  ray-tracing  to  develop  a  phase-error  function  for 
each  atmospheric  realization.  Statistical  mean  and  uncertainty  values  are  obtained  by  running 
multiple  realizations  for  each  engagement  scenario.  For  small  deviations  from  a  nominal 
geometry,  the  error  variance  for  focal  anisoplanatic  effects  will  scale  as 

°^focal  anisoplanaiism  **  0D/H$o)^  ,  (3-1) 

where  D  and  H  represent  the  aperture  diameter  and  beacon  altitude  respectively.  Reasonably 
accurate  Strehl  estimates  can  be  obtained  from  the  error  variance  results  through  the  use  of  the 
extended  Marechal  approximation. 

A  summary  of  some  recent  numerical  results  for  the  SLC-Day  atmospheric  model  and  the 
SWAT  60-cm  aperture  is  given  in  Figure  3-1.  The  lowest  curve  shows  the  expected  Strehl  as  a 
function  of  beacon  altitude  for  a  single  beacon  geometry  and  illustrates  the  deleterious  effects  of 
focal  anisoplanatism.  As  the  source  altitude  is  reduced,  errors  are  incurred  as  a  result  of  the 
unsampled  turbulence  above  the  beacon  and  the  incorrectly  sampled  atmosphere  below  the 
beacon.  The  use  of  multiple  beacons  improves  the  low-altitude  sampling  and  can,  as  indicated 
by  the  middle  curve,  produce  a  significant  Strehl  improvement.  One  of  the  principal  tasks  of  the 
SWAT  field  experiments  is  to  measure  the  performance  difference  between  the  single-beacon  and 
four-beacon  sampling  geometries. 
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Figure  3-1.  SWAT  performance  curves  for  the  SLC-Day  model. 


Numerical  simulations  of  the  multiple-beacon  reconstruction  process  have  consistently 
demonstrated  that  the  theoretical  performance  limit  cannot  be  achieved  with  a  single-altitude 
beacon  configuration.  Errors  due  to  the  improper  measurement  of  high-altitude,  low-spatial 
frequency  turbulence  components  become  increasingly  more  important  as  the  size  of  the  beacon 
array  is  enlarged.  The  upper-most  curve  in  Figure  3-1  describes  the  performance  of  a  source 
geometry  known  as  hybrid  sampling,  in  which  residual  section-tilts  are  corrected  using  a  low- 
spatial  frequency  measurement  of  the  resonant  return  from  the  earth’s  sodium  layer  at  90  km. 
For  the  SWAT  system,  hybrid  sampling  is  expected  to  yield  a  small  but  measurable  performance 
advantage. 

As  described  in  Section  2,  the  laboratory  tests  were  specifically  designed  to  permit  Strehl 
comparisons  to  be  made  for  a  variety  of  source  arrangements.  Although  the  performance  data 
generated  in  this  manner  may  be  compared  with  the  theoretical  results  just  described,  a  more 
direct  measure  of  the  fidelity  of  the  experiment  is  obtained  by  analyzing  the  specific  atmospheric 
realization  achieved  on  the  optical  bench.  To  this  end,  each  of  the  phase  screens  used  in  the 
experiment  was  mapped  with  the  aid  of  a  Zygo  interferometer,  and  the  resulting  phase  arrays 
were  applied  to  the  numerical  stitching  model. 

Figure  3-2  shows  a  phase  measurement  of  one  of  the  phase  screens  (Screen  #1)  and  its 
associated  anisoplanatic  error  when  placed  at  the  high-altitude  turbulence  position  on  the 
laboratory  bench.  Far-field  beam  profiles  computed  through  a  Fourier  transform  of  wavefront 
functions  for  a  diffraction-limited  beam,  a  point- source  10-km  beacon,  and  an  extended-beacon 
geometry  are  given  in  Figures  3-3a,  b,  and  c  respectively.  The  Strehls  derived  from  the  last  two 
images  are  0.36  and  0.42.  It  will  be  shown  in  the  next  section  that  these  numbers  are  in  excellent 
agreement  with  the  Strehl  measurements  obtained  from  the  laboratory  far-field  camera  data. 
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TILT-REMOVED  FOCAL  ANISOPLANAT1C  ERROR 
Variance  =  1.39  radJ 


Figure  3.3a.  Diffraction-limited  far-field  beam  profile. 


Figure  3  3  b.  Far -field  beam  profile  for  a  single  source  at  10  km 
and  propagation  through  Screen  #1  (Strehl  =  036). 


Figure  3.3  c.  Far -field  beam  profile  using  a  four -source-over age  reconstructor 
with  10-km  beacons  and  propagation  through  Screen  # 1  (Strehl  -  0.42). 


4.  CW  PHASE  COMPENSATION  RESULTS 


High  quality  phase-compensation  data  had  been  generated  on  a  semi-routine  basis  in  the 
SWAT  laboratory  since  the  early  fall  of  1987,  but  the  most  reliable  multiple-beacon  results  were 
obtained  in  the  first  month  of  1988.  During  that  period,  a  concerted  effort  was  made  to  follow  a 
procedure  that  would  allow  direct  comparisons  to  be  made  with  the  computer  predictions 
described  earlier.  To  this  end,  a  laboratory  methodology  incorporating  the  following  measures 
was  instituted. 

•  The  beam  path  was  carefully  shielded  to  reduce  dynamic  turbulence  that  would  be 
uncorrectable  at  the  0.5-2  Hz  update  rates  used  in  these  tests.  In  previous 
measurements  it  was  found  that  this  precaution  was  effective  in  reducing  dynamic 
turbulence  to  a  tolerable  level. 

•  Far-field  background  frames  were  recorded  frequently. 

•  Phase  screen  transmission  measurements  were  made  to  allow  for  accurate  throughput 
normalization. 

•  The  effective  system  focal  length  into  the  far-field  camera  was  computed  by 
measuring  the  spot  motion  as  a  function  of  input  wavefront  tilt. 

•  Frequent  closed-loop  measurements  were  made  without  phase  screens  in  the  beam 
path  to  establish  a  static-aberration  calibration. 


All  phase  screens  were  masked  to  ensure  the  use  of  the  same  sections  throughout  the 
tests,  and  the  screens  were  independently  characterized  with  a  Zygo  interferometer. 


Most  of  the  tests  performed  involved  the  placement  of  Screen  #1  at  either  the  low-altitude  or 
high-altitude  turbulence  position.  The  matrix  of  experimental  parameters  is  summarized  in  the 
table  below. 


TABLE  4-1 

Experimental  Parameters 


Screen  #2 
Position 

CCD  Camera 
Source 

Phase  Sensor 
Source 

Reconstruction 

Algorithm 

Not  Placed 

Collimated 

Collimated 

Single  Source  (conventional) 

Low  Altitude 

1  at  10  km 

4  Sources,  Tilt  Included 

High  Altitude 

4  at  10  km 

4  Sources,  Tilt  Removed 

4-Source  Average 

The  last  column,  containing  the  reconstructor  algorithm  listing,  is  actually  incomplete,  in  that  it 
does  not  explicitly  include  all  of  the  variations  tested  during  this  period.  It  was  determined  early 
in  these  tests  that  reconstructors  that  ignored  section-tilt  errors  performed  poorly,  as  did  those 
with  small  section-overlap  regions;  performance  with  these  obviously  inferior  reconstructors  will 
not  be  reported. 


Before  proceeding  with  the  final  data  presentation,  it  is  of  interest  to  describe  some  of  the 
factors  that  would  tend  to  limit  the  achievable  Strehi.  Some  of  the  more  obvious  sources  of  non- 
correctable  error  are  listed  below. 

•  poor  optical  alignment 

•  uncorrectable  common-path  static  aberrations 

•  non-common -path  phase  aberrations 

•  insufficient  signal  into  the  wavefront  sensor 

•  inadequate  mirror  update  rate 
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Considerable  care  was  taken  in  the  optical  alignments,  which  were  performed  daily,  and  the  use 
of  optical  intensifiers  insured  a  high  signal  level  into  the  wavefront  sensor.  The  mirror  update 
rate  was  admittedly  low  (0.5-2  Hz),  but  the  optical  path  was  well  shielded  from  external 
disturbances.  Of  perhaps  greater  concern  is  the  strong  static  astigmatic  error,  indicated  in  Figure 
4-1,  that  was  found  in  the  common-path  optics,  and  the  uncertain  quality  of  the  non-common 
path  optics  into  the  CCD  camera.  The  cumulative  effect  of  these  errors  is  difficult  to  quantify  at 
this  time. 


Detailed  analysis  was  performed  on  a  set  of  12  data  sessions,  each  containing  at  least  30 
mirror  update  experiments.  Strehl  ratios  were  computed  in  accordance  with  the  techniques 
described  in  the  Appendix  to  this  report.  For  the  optical  arrangement  in  which  the  phase  sensor 
viewed  the  collimated  source,  the  far-field  magnification  constant  was  determined  to  be  3.35 
pixels  of  displacement  per  wave  of  tilt  over  the  full  aperture.  This  factor,  applied  to  Equation 
(A-9),  yields  the  following  Strehl  expression 

Z  (peak )-  Zfe (peak) 

S=-  14.3——  -  —  :  collimated  source  ,  (4-1) 

y)  ~Zb(x,y)J 


where  Z  and  Z^  represent  the  signal  and  background  far-field  image  data.  The  magnification 
constant  was  slightly  modified  when  the  main  focus  element  was  adjusted  to  allow  the  phase 
sensor  to  view  the  10-km  sources.  Since  this  action  changed  the  focus  in  the  common  path,  a 
compensating  adjustment  at  the  CCD  camera  was  also  necessary.  Although  a  direct  measurement 
of  tilt  sensitivity  was  not  made  subsequent  to  these  adjustments,  an  estimate  based  on  a 
comparison  of  well-compensated  beam  profiles  indicates  that  the  effective  focal  length  into  the 
CCD  camera  was  increased  by  about  19%  ,  so  that 


S  *  20.1 


Z(peak )  -  Zfe(peak) 
yy  [z  (x,  y)  -  Zb(x,  y)] 


:  10  km  source  . 


(4-2) 
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The  performance  limit  for  the  system  can  be  determined  from  the  far-field  data  collected 
with  the  non-aberrated  optical  path  and  the  collimated- beam  reference.  The  measured  Strehls  for 
this  arrangement  are  plotted  as  a  function  of  mirror  update  iteration  in  Figure  4-2.  (The  time 
interval  between  experiments  was  approximately  0.75  sec  and  a  loop  gain  of  0.25  was 
employed.)  The  steady-state  Strehl  is  seen  to  reach  a  value  of  just  under  0.8  near  the  end  of  the 
sequence,  which  corresponds  to  an  wavefront  error  of  about  Xj  13.  During  this  same  period,  the 
tilt-removed  static  mirror  figure  was  found  to  be  constant  to  better  than  V40.  These  observations 
are  consistent  with  either  a  large  non-common  path  error  or  the  system's  inability  to  compensate 
for  the  static  aberration  shown  in  Figure  4-1 . 

A  performance  summary  for  the  closed-loop  laboratory  experiments  is  presented  in  Table 
4-2.  For  each  of  the  listed  configurations,  a  Strehl  prediction  taken  from  Table  2-2  is  given  in 
the  seventh  column,  and  a  measured  Strehl  representing  an  ensemble  average  over  at  least  15 
experiments  is  shown  in  the  last  column. 

The  first  three  rows  of  Table  4-2  relate  to  configurations  in  which  the  phase  sensor  viewed 
the  collimated  beam  while  the  position  of  phase  Screen  #1  was  varied.  Focal  anisoplanatic 
effects  will  not  be  observed  in  these  cases,  and  only  losses  due  to  fitting  error  and  scintillation 
should  be  important.  An  excess  phase  error  of  the  order  of  V13  is  observed  for  all  of  the 
collimated  beacon  experiments. 
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TABLE  4-2 

Summary  of  Closed-Loop  Results 


File 

Configuration 

Screen 

Errors 

Ac2 

Cum. 

a2 

1654 

Star  Source 

None 

None 

0.00 

0.00 

1904 

Single-Source 

Reconstructor 

Low 

Fitting 

0.08 

0.08 

1801 

High 

<£ 

0.08 

0.16 

1709 

Single  Source 

None 

None 

0.00 

0.00 

1917 

(10  km) 

Single-Source 
Recons  true  tar 

Low 

Fining 

0.08 

0.08 

1813 

High 

Focal 

0.08 

1.50 

1.66 

1725 

4  sources 

None 

None 

0.00 

0.00 

1928 

(10  km) 

Low 

Fitting 

0.08 

0.08 

1830 

Average  of 
Gradients 

High 

Focal 

0.08 

1.25 

1.41 

1720 

4  sources 

None 

None 

0.00 

0.00 

1924 

(10  km) 

Low 

Fitting 

0.08 

0.08 

1824 

Tilt-Removed 

Stitching 

High 

_2 

°X 

Focal 

0.08 

1.15 

1.31 

Strehl 


Theory  Measurement 


1.00 


0.92 


0.85 


1.00 


0.92 


Measurement  error  is  estimated  to  be  ±5% 

++ 

Simulation  result  using  measured  phase  error  for  Screen  #1 
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The  last  nine  rows  of  Table  4-2  provide  an  overview  of  the  system's  operation  when 
simulated  Rayleigh  beacons  were  positioned  at  a  scaled  altitude  of  10  km.  Six  of  these 
configurations  were  either  free  of  intentional  aberrations  or  employed  a  single  phase  screen  near 
the  pupil  plane.  Focal  anisoplanatic  effects  should  be  negligible  under  these  conditions,  and, 
indeed,  the  measured  Strehls  are  found  to  be  essentially  identical  to  those  obtained  when  the 
wavefront  sensor  viewed  the  collimated  reference  laser. 

The  more  interesting  comparisons  of  performance  achieved  in  the  presence  of  focal- 
anisoplanatic  distortion  induced  by  high-altitude  turbulence  are  addressed  in  the  summaries  for 
data  files  SWAT01 161813,  1830,  and  1824.  As  discussed  in  the  last  section,  Strehl  predictions 
for  the  first  two  of  these  configurations  were  obtained  from  numerical  ray-trace  simulations 
incorporating  the  same  beacon  geometries,  screen  placements,  and  reconstruction  algorithms 
used  in  the  experiment.  In  addition,  the  phase  distortion  used  in  the  program  was  derived  from  a 
Zygo  measurement  of  Screen  #1.  The  agreement  between  the  computer  results  and  the 
laboratory  measurements  appears  to  be  excellent.  It  should  be  noted,  however,  that  the 
variations  in  Strehl  over  the  sampled  ensemble  were  large  (see  Figure  4-3),  corresponding  to  an 
uncertainty  of  approximately  ±0.05. 

Although  the  data  listed  in  this  table  were  generated  from  just  a  single  turbulence 
realization,  both  the  strength  and  spectral  characteristics  of  this  phase  screen  are  thought  to  be 
representative  of  a  typical  atmosphere.  All  of  the  low-altitude  beacon  results  are  judged  to  be 
consistent  with  earlier  predictions  based  on  the  SLC-day  turbulence  model,  and  the  difference 
between  the  single-beacon  and  four-beacon  Strehls,  although  modest,  was  not  unexpected.  In 
view  of  these  results,  it  is  likely  that  much  of  the  work  to  be  conducted  at  Maui  will  be  devoted  to 
Rayleigh  scattering  from  altitudes  below  5  km,  where  the  effects  of  focal  anisoplanatism  are 
more  significant. 
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Figure  4-3.  Strehl  plots  for  Sessions  SWATOI 161813  and  SWAT01 161824  (loop  gain  =  0.5). 


5.  PULSED  COMPENSATION  PERFORMANCE 


The  beacon  lasers  to  be  used  in  the  SWAT  II  tests  will  operate  at  a  repetition  rate  that  is 
much  lower  than  the  atmospheric  time  constant,  thus  making  it  impossible  to  perform  closed- 
loop  phase  correction.  The  returns  from  each  pulse  must  be  accurately  sensed,  manipulated,  and 
transformed  into  phases  that  can  be  applied  to  the  deformable  mirror.  In  this  mode,  calibration 
errors  become  much  more  important  and  may  even  become  the  dominant  source  of  error.  For 
this  reason,  the  SWAT  equipment  was  designed  to  be  easily  and  frequently  calibrated  and 
provision  was  made  to  accommodate  non-linear  conversion  functions. 

In  order  to  test  the  functionality  of  the  system  in  its  pulsed  or  "go-to"  mode,  the 
deformable-mirror  loop  gain  was  set  to  unity,  and  far-field  camera  measurements  were  made 
immediately  after  a  single  mirror-updare  cycle.  The  timing  sequence  was  similar  to  that  to  be 
used  in  the  field,  allowing  only  a  millisecond  for  the  completion  of  the  phase  sampling, 
reconstruction,  and  mirror  correction.  A  performance  summary  comparing  the  pulsed  and 
closed-loop  modes  is  shown  the  table  below  and  in  Figure  5-1.  All  of  the  data  were  generated 
using  a  single  collimated  reference  beacon. 
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TABLE  5-1 


Pulsed  Operation  Summary 
(Files  1949,  1952,  2003,  &  2005) 


Closed-Loop  Operation 

Pulsed  Operation 

Turbulence 

Conditions 

rms  Phase 
Fluctuation 
(waves) 

Strehl 

rms  Phase 
Error 
(waves) 

Strehl 

Static  Aberration 

0.04 

0.65 

0.08 

0.63 

Phase  Screen  #1 

0.05 

0.65 

0.08 

0.60 

A  comparison  of  the  closed-loop  and  pulsed-mode  results  shows  the  performance 
estimates  to  be  nearly  identical.  Since  these  test  were  conducted  under  relatively  strong 
turbulence  conditions,  they  represent  a  fairly  convincing  end-to-end  demonstration  of  pulsed- 
mode  functionality.  This. experiment  will  undoubtedly  be  repeated  during  the  field  integration 
and  testing  process. 
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6.  SUMMARY 


The  Lincoln  Laboratory  SWAT  program  represents  a  unusually  ambitious  attempt  to  gauge 
the  effectiveness  of  a  diverse  set  of  synthetic-beacon  concepts.  The  laboratory  tests  conducted  in 
January  1988  were  intended  to  exercise  the  key  elements  of  the  adaptive-optics  subsystem,  the 
data-collection  equipment,  and  the  data-analysis  software.  It  was  also  hoped  that  these 
experiments  would  provide  an  initial  performance  evaluation  of  single  and  multiple-beacon 
phase-compensadon  techniques  that  could  be  quantitatively  compared  with  current  theory. 

In  all  of  these  aspects  the  experiments  are  judged  to  have  been  successful.  In  particular,  the 
agreement  between  theory  and  experiment  in  the  closed-loop  operational  mode  was  excellent  in 
essentially  every  beacon  configuration.  Successful  correction  was  obtained  in  the  single-pulse, 
"go-to"  mode,  with  resulting  Strehl  ratios  only  slightly  worse  than  obtained  with  CW,  closed- 
loop  compensation. 
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APPENDIX: 

FAR-FIELD  IMAGE  ANALYSIS 


The  processing  methods  used  to  extract  estimates  of  Strehl  from  far-field  imagery  are 
typically  reviewed  at  the  beginning  of  each  new  field  exercise,  and  the  SWAT  experiments  have 
been  no  exception.  Over  the  years  a  wide  variety  of  approaches  have  been  developed  to  deal 
with  this  problem,  all  of  which  display  a  fairly  high  susceptablity  to  data  normalization  errors. 
Following  a  careful  preliminary  analysis  of  the  data  obtained  during  the  SWAT  laboratory 
exercise,  a  new  approach  has  been  devised  that  appears  to  be  somewhat  better  suited  to  this  task. 
The  process  comprises  two  distinct  stages:  during  the  First  step  the  raw  data  are  carefully 
corrected  to  remove  fixed-pattern  noise  and  stray  light  effects;  in  the  second  step  an  estimate  of 
the  peak  irradiance  value  is  obtained  from  an  energy  distribution  function.  The  essential  elements 
of  this  process  are  outlined  below. 

The  high-reso^tion  CCD  cameras  used  in  the  SWAT  program  can  provide  two- 
dimensional  beam  images  of  extremely  high  quality.  The  raw  output  data  Z(x,y)  can  be  modeled 
as  follows 

Z(x,y)  =  xd  x0  d2  g(x,y)  E(x,y)  +  Zb(x,y)  ,  (A-l) 

where  E  and  Z  represent  the  input  optical  irradiance  and  the  electrical  output  signals  respectively, 
is  the  camera  dwell  time,  tq  is  the  optical  throughput,  d  is  the  pixel  dimension,  g  is  the  the 
camera  gain,  and  Zb  is  the  fixed-pattern  background.  The  background  is  normally  assumed  to  be 
spatially  varying  but  temporally  static  so  that  the  Zb  term  can  be  eliminated  by  subtracting  camera 
frames  taken  with  the  input  light  blocked.  It  is  further  assumed  that  the  gain  is  independent  of 
pixel  location  but  that  it,  and  the  system  throughput,  cannot  be  accurately  measured.  Therefore, 
the  optical  irradiance  function  at  the  camera  focal  plane  bears  the  following  relationship  to  the 
output  data 

E(x,y)  =  (1/y)  [Z(x,y)  -  Zb(x,y)]  (A-2) 

where  y  is  indeterminate.  Therefore,  a  direct  radiometric  determination  of  the  Strehl  ratio  is 
precluded  but  a  reasonably  accurate  estimate  can  be  obtained  if  the  diffraction-limited  beam 


A-l 


profile  for  the  collection  aperture  can  be  computed.  In  general,  a  parameter  having  the  form 

Z  (x,  y )  -  Z  bU,y) 

££[ZU,y)  -  Z  (x.y)  j 

x  y 

will  be  generated  at  some  point  in  the  process,  which  eliminates  the  y  term.  This  energy 
normalizauon  is  less  straightforward  than  it  appears,  owing  to  the  difficulty  of  simultaneously 
achieving  high  spatial  resolution  across  the  beam  while  obtaining  an  accurate  measure  of  the  total 
flux.  Separate  narrow  and  wide  field-of-view  detectors  can  be  used  for  this  purpose,  but  care 
must  be  taken  to  establish  an  accurate  mutual  calibration  of  the  two  sensors. 

A.l  BASELINE  CORRECTION 

Failure  to  properly  subtract  the  fixed-pattern  background  noise  is  the  greatest  potential 
source  of  error  in  the  Strehl  calculation,  since  a  baseline  offset  on  the  order  of  0.1%  of  the  peak 
pixel  intensity-  can  translate  into  a  10%  error  in  the  estimate  of  the  total  beam  energy.  An  offset 
error  will  occur  when  the  ambient  lighting  changes  between  the  data  and  background  frames  or 
when  the  camera  is  exposed  to  signal  light  scattered  from  optical  elements  within  its  field  of 
view. 

Recall  that  the  Z  -  Zb  operation  is  intended  to  produce  a  data  field  that  is  free  from  fixed- 
pattern  errors.  Assuming  that  the  focused  beam  subtends  only  a  small  fraction  of  the  total 
detector  array  area,  a  histogram  of  Z  -  Z*,  should  display  a  Gaussian  peak  centered  at  zero.  Any 
displacement  of  this  peak  from  the  null  position,  as  illustrated  in  Figure  A-l,  is  evidence  of  a 
baseline  offset  that  should  be  corrected. 

Baseline  compensation  emails  the  formation  of  the  quantity  l(x,y)  =  Z(x,y)  -  Zb(x,y)  -  Z0, 
where  Zq  is  the  computed  mean  of  the  baseline  peak  of  the  Z  -  Zb  irradiance  histogram.  While 
the  absolute  value  of  Zq  is  typically  small,  the  Strehl  measurement  error  is  proportional  to  the 
product  of  Zr,,  and  the  number  of  pixels  in  the  detector  array.  The  implementation  of  this  simple 
proceuure  has  significantly  improved  the  frame -to -frame  uniformity  of  the  Strehl  computations 
performed  on  the  SWAT  laboratory  data. 


FREQUENCY  of  OCCURRENCE 


A.2  PEAK  IRRADIANCE  ESTIMATION 

For  a  uniformly  illuminated  circular  aperture,  the  diffraction-limited  beam  profile  is 
specified  by  the  airy  irradiance  function 


( (kD/2f)2' 

f  J.(kDr/2f)l 

v  4lt  7 

L“  kDr/2f  J 

(A-3) 


where  D  is  the  aperture  diameter,  f  is  the  focal  length  of  the  telescope,  and  Q  is  the  pulse  energy. 
It  will  be  convenient  to  approximate  this  function  by  the  Gaussian 
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where  Oq 


defines  the  Gaussian  beam  width. 


The  beam  profile  of  an  aberrated  beam  might  be  modeled  in  the  following  way 


E(r)  = 


Qs 
2k  o' 


exp 


Ac2) 


(A-5) 


where  S  is  the  peak  irradiance  (the  usual  definition  of  Strehl)  and  c  is  the  beam  diameter  (which 
will  usually  be  larger  than  aQ).  The  irradiance  can  be  treated  as  a  probability  density  from  which 
a  distribution  function  can  be  formed 


F(r)  =  2k J  E(r)  rdr  =  Q  S 
0 


1  -  exp 


11  r 

2 


2  3 


a2,j 


(A-6) 


and 


F(E)  =  QS 


r^\-2KC?E 

VaoJ 


(A-7) 


The  function  F(E)  represents  the  total  focal  plane  energy  attributed  to  those  pixels  exceeding  the 


irradiance  threshold  E.  F(E)  is  seen  to  be  a  linear  function  of  the  irradiance,  with  a  slope  that  is 
proportional  to  the  beam  radius  0.  The  x-axis  intercept  Ep  =  QS/2naQ2  provides  a  means  to 
measure  of  the  peak  irradiance,  while  the  y-axis  intercept  gives  the  total  energy  included  in  the 
Gaussian  core.  Therefore,  a  linear  fit  to  the  upper  portion  of  the  irradiance  distribution  function 
will  generate  estimates  of  both  the  Strehl  and  the  effective  beam  diameter.  The  general 
appearance  of  the  F(E)  curves  for  diffraction  limited  and  aberrated  beams  are  illustrated  in  Figure 
A-2.  The  knee  of  the  curve  will  generally  occur  well  below  the  Ej/2  point,  so  it  is  usually  safe  to 
include  all  data  points  above  this  value  in  the  linear  regression  operation. 


The  preceding  derivation  provides  a  working  definition  of  the  Strehl  ratio  that  includes  the 
following  analysis  steps: 

1 .  the  formation  of  the  corrected  camera-data  histogram  h(f), 

2.  the  creation  of  the  camera  distribution  function  F(I)  =  ^  i  h(i) , 

i>I 

3.  the  development  of  a  linear  fit  to  the  upper  half  of  the  F(I)  curve  to  obtain  F(I),  and 

4 .  the  determination  of  the  x-axis  intercept  Ip  for  which  F(lp)  =  0. 

The  Strehl  ratio  is  now  defined  to  be 

2*0^ 

S=IIl(x,y)  •  <A-8> 

*  y 

As  usual,  the  determination  of  Strehl  requires  an  accurate  measurement  of  the  total  pulse  energy 
and  knowledge  of  the  camera’s  magnification.  Observe  that 


(A-9) 


where  d  is  the  pixel  dimension,  n  is  the  number  of  phase  sensor  subapertures,  and  A  is  a 
magnification  factor  given  in  units  of  beam  motion  (pixels)  per  wave  of  tilt  per  subaperture. 
Since  the  Strehl  is  proportional  to  A2,  it  is  important  that  this  magnification  factor  is  remeasured 
following  each  change  to  the  optical  arrangement. 


A-5 
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Figure  A2  Energy  distribution  functions  for  diffraction-limited  and  aberrated  beams. 
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